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ON | The chiral surface electrons in the bulk quantum Hall effect probably form the first extended 

system in which conductance fluctuations can be calculated non-perturbatively in the presence of 
disorder. By use of the Kubo formula with appropriate boundary conditions, we calculate exactly 
the variance of conductance with non-perturbative methods. We find that the conductance fluctu- 
ations of this system are nonuniversal and the variance of the conductance scales in a very peculiar 
way. This result can be checked with exact computation using the Landauer-Buttiker formula and 

on : both methods show the same scaling behavior. We have also calculated the diffusion constant fluc- 

tuations exactly. We find that the diffusion constant fluctuations vanish and thus play no role in 

| | the conductance fluctuations. 

43 : 

Electronic systems under the influence of disorder exhibit many intriguing phases most of which are still not 
well understood [Q . The understanding of the transition of electronic states from localized to extended is still 
far from complete. In the standard scaling theory of localization j3|, it is argued plausibly that the beta function 
'"" J ■ P(g) = dlng/d\nL is a function of the dimensionless conductance g alone. This one-parameter scaling idea received 
further support in the metallic regime from the experimental appearance of universal conductance fluctuations (UCF) 
[0] in small metals as well as from an interesting numerical simulation |^] which revealed the UCF as a consequence 
of quantum interference. When the electronic states are localized, however, the system generally has a very broad 
conductance distribution |l| . The average conductance thus might not seem enough to fully characterize the system. 
Indeed, many-parameter scaling Q has been proposed to take into account the broad conductance distribution 
near the derealization transition. There is no method to determine from first principles whether one-parameter 
or many-parameter scaling is more appropriate except by calculating the higher cumulants of the conductance, such 
as conductance fluctuations (CF), directly. So far CF calculations have been done only perturbatively |t],|| and 
are valid only in the metallic regime. It is thus natural to ask will these perturbative results hold to all orders of 
. perturbation and what happens in the regime that is inaccessible to perturbative methods. An exact conductance 
£C) ' calculation in any physical system is thus important. 

We now introduce an interesting 2D anisotropic system, the chiral surface states of the bulk (multi-layered) integral 
quantum Hall effcct(IQHE). In the multi-layered quantum Hall sample, at the edge of each layer the electrons circulate 
in one direction only and can therefore be modeled as chiral Fermions J9| . When the tunneling of these edge electrons 
in between the layers is allowed, the edge states then smear over a two dimensional sheet where electrons move 
ballistically transverse to the field and diffusively otherwise. Bechgaard salts and multi-layer hetrostructures Jl|] 
are candidates for this system. The collection of these edge states turned out to be a novel metallic phase and thus 
is also termed a chiral metal J13| . This system is critical in the sense that the only length scale of the system is the 
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system size 12 . One peculiar feature of the chiral metal is that it can retain conductance much smaller than e 2 /h 
while still scaling ohmically. This should be contrasted with the ordinary metal for which ohmic scaling only happens 
q ' at large conductance. 

We shall focus on the limit in which the length along the chiral direction (circumference) L x is much greater than 
the transverse dimension L z ; this corresponds to the zero-dimensional regime defined by Balents et al. The 
finite L x case will be discussed in a future publication ]l5| . Starting with linear response theory, we derive the Kubo 
formula suitable for this system and calculate the disorder-averaged conductance as well as perform the first exact 
conductance fluctuation calculation non-perturbatively. We find that the CF are non-universal because the variance of 
the conductance is a function of the hopping amplitude t, which in turn is related to the system's diffusion constant. 
The CF scale with the system size as L X /L 2 Z . We also find that the disorder-averaged conductance scales ohmically: 
(g) cx L x /L z where () is used to denote the disorder average. Define 5g = g — (g) and A = (5g 2 )/(g) 2 . We have 
A ~ l/L x <C 1 and thus g is narrowly distributed. Furthermore, A — -» as L x — ► co means that g is self-averaging. 
We have also calculated the average conductance using the Landauer-Buttiker(LB) formula as well as the Einstein 
relation. These two results are identical. We have also calculated the diffusion constant fluctuations exactly and find 
surprisingly that there are no diffusion constant fluctuations at all. This should be contrasted with the conventional 
result for ordinary metals with dimensionality two or lower for which the major contribution of CF is the diffusion 
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fluctuation! This can be viewed as another salient feature of the novel chiral metallic phase. We can also calculate the 
CF by using LB formula, but the sums involved are difficult and are thus done numerically. The CF obtained from 
the LB formula agree qualitatively with the Kubo calculation and show the same scaling behavior. One interesting 
feature worth mentioning is that the CF vanish when t = or t — 2 which correspond to the situations of insulating 
and maximum-conducting respectively. We find this happens both in LB type of calculation as well as in the Kubo 
type of calculation. 

For simplicity, we consider the situation where only the first Landau level is filled, i.e. only one edge state per 
quantum Hall layer. The system thus consists of a collection of 1 + 1 dimensional ballistic states moving in the positive 
x (chiral) direction, each described by a chiral Fermion, which are coupled together by a nearest-neighbor hopping 
amplitude t(n,x): 



H =^2j dx -ipl id x ip n - [t(n,x)iplil) n+1 + h.c.J 



(1) 



There are two versions of the chiral model distinguished by the form of the hopping amplitude t(n,x). The first one 
we shall call continuous hopping has t(n, x) = t being a constant. The second one which we shall call discrete hopping 
has t(n, x) as a sharply peaked function at even(odd) integer x for n being even(odd) and t(n, x) being zero otherwise. 
We shall focus only on the discrete hopping case. The chiral fermion velocity v is 1 in our choice of units. The full 
Hamiltonian including random scattering is Ti. = Tio + where 

H!=J2 [ dxV n (x)i>ti> n . (2) 

n J 

An illustration of this system is given in Fig. 1. The random potential V n (x) is assumed to have zero mean. To make 
our discrete hopping model identical to the one used in ref. , we require that the integral of the random potential 
over each link J,. , V n {x)dx mod (2-7r) is a random variable distributed uniformly over (0, 2tt\. In the large L x limit, 
due to the ballistic nature of the electron motion along the x (chiral) direction the Green's function which is retarded 
in time must also be retarded in x The retarded Green's function with constant frequency E, R E , satisfies the 
following equation 



id x + V2n(x) — E j R E (2n, x; m, x') + k R E (2n — 1, x; m, x') + k e R E (2n + 1, x; m, x') = — 62 n ,mS( x ~ x> ) 

( - id x + V 2n -i(x) - E^R E (2n - l,x;m,x') + h R E (2n, x; m, x') + A e R E (2n - 2,x;m,x') = -S 2n -i,m5(x ~ x') (3) 

where A Q ( e ) (x) = t y^° dd ( cvcn ) j^ x _[^ a nd f(x) is a symmetric function sharply peaked around zero with the condition 
that / f(x)dx — 1 p7[ . Note that the nonzero fixed frequency E can be gauged away by a simple local gauge 
transformation G — > e X G, 

The retarded Green's function R(n, x; n', x') obeying Eq. (^) can be expressed as a sum over paths of electron going 
from (n' ,x') to (n,x) 

R{n, x- n',x') = J2l[ [xi W exp(i0j)] . (4) 
p lep 

As depicted in Fig. 1, the electron only moves upwards along the links and can hop horizontally to its neighboring 
layer across the nodes. Each path contributes an amplitude which is a product of the phases (f> — J link V n (x')dx' of 
the constituent links together with a prefactor \(t) for each node which depends on whether the fermion hops to a 
neighboring layer or remains in the same layer. The prefactors can be obtained by integrating Eq. (^) across a node 
. An analogous expression can be written for the advanced Green's functions. 
Let us first consider the disorder-averaged conductance which will illustrate the methods used here to calculate the 
CF. The disorder-averaged conductance has previously been calculated for this model by Chalker and Dohmen Jl2| . 
Using the standard linear response theory a Kubo formula can then be written that expresses the conductivity in 
terms of the Green's function as 



cr zz (n, x; n' , x') = — A n _i(x)A„/_i(a;') I [R(n — 1, x; n' — 1, x')A{n' , x'; n, x) + R(n, x; n' , x')A{n' — l,x'; n — 1, x) 

- R(n, x; n' - 1, x')A(n' , x';n-l,x)- R(n -l,x; n' , x')A{n' -l,x'; n, x)] + [R <->• A] |e 2 t 2 /47r (5) 
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where the E = subscript is suppressed from the Green's functions. In deriving the above expression it is important 
to keep in mind that the current operator along the non-chiral direction is J n (x) = ietA„__i(a;)[^^ n _i — ^l-iV-Vi] 
where A n (x) = A e (x)/t if n is even and A n (x) — A (x)/t if n is odd [ p^[ . 

Using the path integral expressions for the retarded and advanced Green's function, Eq. (||) can be expressed as 
a double sum over paths of retarded Green's function multiplied by the paths of advanced Green's function. Upon 
disorder-averaging, due to the random link phases, only the diagonal terms, in which the retarded and advanced paths 
pair up, will survive in this double sum. Hence, the disorder-averaged conductivity is nonvanishing only when both 
x and x' are equal to the same integer and n — n' Jl8| 

(a zz (n, x] n', x 1 )) = e 2 t 2 /2ir 8 nt7l >8 x , x >. (6) 

This leads to (g) = (L x /2L z )e 2 t 2 /2ir; the factor L x /2 comes from the fact that A e ( Q ) only sum over even(odd) integers 
along the chiral direction. 

The disorder-averaged conductance can also be evaluated using the LB formula. This approach differs from the Kubo 
method in that scattering boundary conditions at the probes are implemented exactly. To proceed we have to calculate 
the probability 2 of finding outgoing current at z = L z , x = j while constant current is injected at z = 1, x = i. The 
conductance is given by g = yV J —1 |ij,j| 2 . Note that in terms of diffuson D(n,x;n' ,x') = R(n, x; n', x')A(n', x'; n, x), 
\tij\ 2 is given by D(L z ,j; 1, i). Calculating (g) is thus equivalent to counting paths with appropriate weights. A direct 
calculation thus shows that 

, x _e 2 t 2 {L x /2) 

W 2tT< 2 + (1-< 2 /4) 2 L 2 [1 

where the factor L x /2 originates from the fact that along the chiral direction there is only one lead out for every two 
lattice spacings. The above result holds for all t and agrees with the result from the Kubo formula for small t. At 
t = 2 where the system becomes ballistic also in the non-chiral direction, the disorder-averaged conductance saturates 
as one might expect from a Landauer conductance calculation. The result from Kubo formula, however, does not 
show this behavior as t — > 2. It will be interesting to understand the discrepancies between Kubo and Landauer 
results. For t ^ 2 and large L z , {g) ~ e 2 t 2 L x /[A'K{\ — t 2 /A) 2 L Z ] which reduces to the result obtained in ref [|l2| with 
+ t 2 /4) identified as the t used in ref (l|. 
To calculate the average conductance in the third way, let us first consider the diffusion constant D which is defined 
as the large L limit of the quantity ^ n n 2 D{n,L) /2L where D(n,x) = D{n,x;n' = Q,x' = 0). We find that JlJ] 
the disorder-averaged diffusion constant is given by (D) = t 2 /2(l — t 2 /4) 2 + 0(1/L). As t — > 2, it seems that (D) 
diverges. A careful analysis by having L large but fixed shows that as t — > 2 one gets (D) ~ L/2. This signifies that 
when t — » 2 the electron motion becomes ballistic in the non-chiral direction also. We also calculate exactly [Q the 
diffusion constant fluctuations and find them vanish for arbitrary L: 

((D - (D)) 2 ) =J2n 2 m 2 [(D(n,L)D(m,L)) - (D(n, L)){D(m, L))]/AL 2 = 0. (8) 

This indicates that the distribution of the diffusion constant is 5(D — Do) independent of the disorder. If we wish 
to use the Einstein relation to obtain conductivity in this case, we only need to calculate the disorder-averaged 
density of states which is 1/2tt. Assuming Ohmic scaling, we then obtain the disorder-averaged conductance (g) = 
e 2 t 2 L x /[4ttL z (1 — t 2 /A) 2 ] for t ^ 2. This result is exactly the same as the result from using the LB formula at large 
L z and agrees with the result from the Kubo formula for small t. 

The calculations of CF as well as the diffusion constant fluctuations involve disorder-averaging the product of four 
Green's functions. But as already mentioned, since the Green's functions must pair up, we can turn this average into 
a two-diffuson problem with contact interaction which comes from the possibility of exchanging paired-partners when 
two diffuson cross spatially. The quantity of central importance is the probability of finding one diffuson at position 
i"i and the other at position Y2 for two diffusons starting at their specified positions. In the large L x limit where the 
loops circled the circumference once or more can be neglected, we can write down an evolution equation along the 
chiral direction for the interacting two-diffuson problem. We then perform a discrete Laplace transform and introduce 
the center of mass momentum as well as the relative momentum of the two diffusons. The results are obtained in 
Fourier space and are then converted back to real space A simpler version of this calculation, although it appears 
in a different context, has also been done |lj|. 

To calculate the CF exactly via the Kubo formula, we first observe that g — (g) is obtained by throwing away 
terms with x — x' in Eq. Note that the integral of A„_i(x) over x gives rise to ^° vcn ( odd ) for n — 1 even(odd). 
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Furthermore, since the contribution to the CF from x > x' and x < x' will be the same due to the symmetry in Eq. (|a), 
we can work out the CF contributed from x > x' and then multiply it by two. Denote D(n, ri) = D(i_,n; l' + , ri) with 
both I and V being integers. After some calculations [Q, in which the zero wave- vector mode along the non-chiral 
direction of diffuson was taken out, we obtain 

(5g 2 ) = (e 4 i 4 /47rL 4 ) (D(n- 1; ri - l)D(n; ri) + D(n - 1; n')D(n; ri - 1) 

;=1 l'=0n(!),n'(l') 

+ (1 - t 2 /4) 2 t 2 /(l + i 2 /4) 4 [D(n - 1; ri - l)D(n; ri - 1) + D(n; ri)D(n - 1; ra')] ) 

= [e 4 t 2 (l-t 2 /4) 2 (l + r) J L x ]/[4 7 r 3 (l + 2r) J L2] when L x » L 2 > 1 (9) 

where T = t 2 (l — i 2 /4) 2 / (l+i 2 /4) 4 . For t ^ 2 and large L z , we find that {Sg 2 } scales with the system size as L x /L 2 . It 
is interesting to note that only in the t — > limit where perturbative calculation being plausible, our exact calculation 
agrees with the result of a perturbative calculation JlGj . For t = 2 where electronic motion becomes ballistic also 
in the non-chiral direction, we find the CF vanish as physically expected. For general t, CF exhibit complicated t 
dependence which is not accessible to perturbative calculations. 

Starting from LB formula (g 2 ) contains disorder-averaged terms of products of two i;.j| 2 . To show how this is 
done, let us consider the case i' > i and j > f for concreteness. The quantity (|ti,j| 2 |ti',j'| 2 ) now can be written as 

(\Uj\ 2 \U'j>\ 2 ) = £ {D(z 1 ,i';l,i)){D(L z ,j;Z2;j')){D(z2,j';z 1 ;i')D(L z ,f;l,i')). (10) 

Zl,Z2 = l 

Following the methods outlined above, one can calculate the disorder-averaged one-diffuson or two-diffuson explicitly. 
The sums over z\, z%,i,i j, j turn out to be quite complicated. We thus perform the sums numerically. We have 
performed the calculations for L z up to 10 and L x up to 6000. To extract the scaling behavior of the CF, we plot 
(L x / {Sg 2 )) 1 / 2 against L z , and find that it is very much a straight line for most values of t in good agreement with 
the results found in using the Kubo formula. The case for t = 0.31 is shown in Fig. 2. We have also check the specific 
cases where t = and t = 2; and we see that the CF vanish in both cases in agreement with the Kubo calculation. 

To conclude, we have performed the CF non-perturbatively by using the path integral expression of the Green's 
function. We find that the CF are non-universal and scale in a peculiar way. It has been speculated [^0| that UCF 
always appear for mesoscopic system with extended quantum states. Our exact calculation shows that at least this 
system is an exception. 
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FIG. 1. The directed chiral network model. Full and dashed lines represent links and nodes respectively. 
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FIG. 2. The abscissa denotes L z and the ordinate marks the corresponding [L x / (Sg 2 )) 1 ^ 2 . The hopping parameter t = 0.31. 



